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Abstract 

'> ' 

\^ ' Consider the stochastic partial differential equation 

(N' 



a 
> 



Ut = Uxx + U'W, 



(N 

qT I where x ^ I = [0, J], W = W{t, x) is 2-parameter white noise, and we 

Q> ■ assume that the initial function u{0, x) is nonnegative and not identi- 

ty . cally 0. We impose Dirichlet boundary conditions on u in the interval 

I. We say that u blows up in finite time, with positive probability, if 

there is a random time T < oo such that 



P ( lim sup u(t, x) = oo ) > 0. 



X 

C^ ' It was known that if 7 < 3/2, then with probability 1, u does not blow 

up in finite time. It was also known that there is a positive probability 
of finite time blow-up for 7 sufficiently large. 

In this paper, we show that if 7 > 3/2, then there is a positive 
probability that u blows up in finite time. 
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1 Introduction 



We consider the heat equation with a nonhnear additive noise term. 
ut = u^x + u^W, t > 0, x G I = [0, J] 



:i-i) 



M(t,0) 



u{t, J) 
uo(x) 







Here, W = W{t, x) is 2-parameter white noise, 7 > 1, and uo{x) is a continu- 
ous nonnegative function on I, vanishing at the endpoints, but not identically 
zero. Suppose that we are working on a probability space {Q,J-',P), and fix 
a point uj & Q. If there exists a random time T = T{uj) < 00 such that 



limsupM(t, x) 



00 



then we say that u blows up in finite time (for the point u). 

For deterministic partial differential equations, there is a large literature 
about blow-up in finite time. See ||buj6q| , |[b'lV185|| , |[b'M8(J|| , |FK92[, and 



LN92|| for example. Suppose that we are dealing with the equation 

dw(t,x) 



dt 
w(0,x) 



Aw(t, x) + g{w{t, x)) 
wo{x). 



One basic ideas is the following. Suppose that g{w) increases faster than 
linearly. If a high peak forms in the solution w{t, x), then the term g{w{t, x)) 
will win out over the term Aw{t, x), the growth of the peak will be governed 
by the ordinary differential equation 

w'{t)=g{w{t)). 

We can solve this equation explicitly, and its solutions often blow up in finite 
time. 

On the other hand, for stochastic partial differential equations (SPDE) 
there are very few papers about finite-time blow-up. Apart from the heat 



equation, the author ||Mue97a|| has studied the wave equation 
d'^Ut(x) 



duo{x) 
dt 

Uo{x) 



Aut{x)+g{ut{x))W{t,x) 

hi{x) 
ho{x). 



t>0,xeR 



where g{u) grows like u{logu)" for some a > 1. For g{u) = m" with a > 1, 
one would guess that solutions would solutions would blow up in finite time. 
But finite time blow-up is not known for any value of a. Similar techniques 



were used in ||Mue93|| to obtain a modulus of continuity for solutions of the 
wave equation with noise in higher dimensions, with correlated Gaussian 
noise instead of white noise. 

There is more precise information about the heat equation with noise. 
Suppose that -u is a solution to (|l.lj ). In |[Mue91|| it was shown that if 7 < 3/2, 



then, with probability 1, u does not blow up in finite time. Krylov [|Kry94|| 
gave another proof of this fact for a more general class of equations. The 
papers | Mue97a | and [ Mue97b ] are also relevant. We refer the reader to Par- 
doux [ pr'ar93 | for this and other questions about parabolic SPDE. Returning 
to the question of blow-up, it was shown in |[MS93|| that there exists 70 > 1 
such that if 7 > 70, then with positive probability, u blows up in finite time. 
The argument in [|MS93|| was not sharp enough to give the best value of 70, 
and the question of whether one could take 70 = 3/2 was left open. The 
main theorem of this paper answers this question in the affirmative. 

Theorem 1 Let u{t,x) satisfy (\1.^ ), and suppose that 7 > 3/2. Then, with 
positive probability, u blows up in finite time. 

Of course. Theorem |l| does not tell us what happens at 7 = 3/2. Sur- 
prisingly, the proof of Theorem |l| uses many of the same ideas as in |[MS93|| , 
although in a sharper form. 

Now we discuss the rigorous meaning of (|1.1|) , following the formalism of 
Walsh ||Wal86| | , chapter 3. Before giving details, we set up some notation. Let 
G{t,x,y) be the fundamental solution of the heat equation on I. If G{t,x) 
is written as a function of 2 variables, we let G{t, x) be the fundamental 
solution of the heat equation on R. In other words 

G{t,x) = 

It is well known that 

Git,x,y)<Git,x-y). 

We regard ( |1 . 1|) as shorthand for the following integral equation. 




u{t,x)= G{t,x,y)uo{y)dy+ G{t - s,x,y)g{u{s,y))W{dyds) (1.2) 

Ji Jo Ji 



where the final term in ( |1.2| ) is a white noise integral in the sense of [ Wal86 



Chapter 2. Because g{u) is locally Lipschitz, standard arguments show that 
( |1.1| ) has a unique solution u{t,x) valid up to the time 0"^ at which |M(t,a;)| 
first reaches the level L for some a; G I. Similar arguments are given in 
| Wal86|] , Theorem 3.2 and Corollary 3.4, and his reasoning easily carries over 



to our case. Letting L —^ oo, we find that ( p..l|) has a unique solution for 
t < a, where a = lim^^oo c"l- H cr < oo, one has 

lim sup \u{t, x)\ = cxD. 

Our goal is to show that a = oo with probability 1. 
More generally, we regard 

vt = Vxx + g{v)W, t>0,xel 
v{t,0) = v{t,J) = 
f(0,x) = vo{x) 

as a shorthand for the following integral equation, which may only be valid 
up to some blow-up time. 

v{t, x)= G{t, X, y)vo{y)dy + G{t - s, x, y)g{v{s, y))W{dyds) 

Ji Jo Ji 

Lastly, we will always work with the cr-fields Tt = T^ generated by the 
white noise up to time t. That is, Tt is the cr-field generated by the random 
variables /q /j 0(s, x)Vr((ix(is), where varies over all continuous functions 
on [0, t\ X I. 

We now summarize the argument in [|MS93|| , which is based on the anal- 



ysis of the formation of high peaks. Such peaks will occur with positive 
probability. We wish to show that, with positive probability, such peaks 
grow until they blow up in finite time. If a high peak forms, we rescale the 
equation and divide the mass of the peak into a collection of peaks of smaller 
mass, and these peaks evolve almost independently. In this way we com- 
pare the evolution of m to a branching process. Large peaks are regarded as 
particles in this branching process. Offspring are peaks which are higher by 
some factor. We show that the expected number of offspring is greater than 
one when 7 > 3/2, and thus the branching process survives with positive 
probability, corresponding to blowup in finite time. 

Finally, we remark that in ( |1 . 1|) , we could replace u' with a function g{%i) 
satisfying giu) > cu"' for some c > 0. Then Theorem ^ would still hold, 
provided 7 > 3/2. 



2 Proof of Theorem [1] 

We give a proof by contradiction. Assume that 

P{a <oo)=0. (2.1) 



We recall Lemma 2.4 of ||MS93|| 



Lemma 1 Suppose that u solves (\1.H) up to some T^ stopping time r. Let 
L> 0. If we let 

t > 0, xe iZ^^^-^) 





v{t, 


x) = L~\ 


(tL^^-^^ 


xZ^d 


-') 


then v{t, x) 


solves 












dv 
dt 


dx^ 


+ b{v, 


i)W 






v{t,0) 


= v(t, 


jr^^- 


^.). 






v{0,-) 


= Vo 


t> 0, ( 





< X < Jl2(7-1) 

up to the J-''^ -stopping time rL'^'^^^'^\ where vq{x) = L^^Uq (xL'^^^^'^^j for all 

< X < JU^^<-'^\ for allt >0 andO < X < jD^'^-'^\ and W is the white 
noise on S(R+ x [0, JL^^^^^^]) defined by 

W{A) = L'^^-^^ I XA{tL^^^-^\xL^^^-^^)W{dt, dx) 

JR+xI 

for all A in 'B(R+ x [0, JL'^'^'^^^^]) with finite Lebesgue measure. 



In | MS93|| we fixed L, which we called L, and took 7 to be very large. In 



the current proof we wish to deal with all 7 > 3/2, so we take L as our large 
parameter. We will find that the probability of a peak getting up to level L 
is about p = 1/L, from the gambler's ruin problem. Using Lemma [1|, we will 
see that after rescaling, a peak of size L gives rise to A^ = L'^^"'~^^ offspring. 
Thus, the expected number of offspring of our initial peak should be 

pN = (1/L) ■ L2(^-i) = L'^~\ 

If 7 > 3/2, then 27 — 3 > and pN —>■ 00 as L —>■ 00. Of course, the above 
heuristic calculation will suffer from the rough estimates we make during 



the course of the proof. Our hope is that taking pN is large enough will 
compensate for all of our sloppiness. 

We will consider solutions u{t, x) to a slightly more general equation than 



o 



ut = u,, + g{u)W, t>0,xeI=[0,J] (2.2) 

u{t,0) = u{t,J)=0 
u{0,x) = Uo{x) 

where g : [0, oo) — > [0, oo) is a locally Lipschitz function satisfying g{0) = 



and g{u) > u^ for u > 0. The same argument as for (|1 . 1|) gives existence and 
uniqueness of u{t, x) up to the blow-up time for u. 

Let (/^(t, x) = (f^'^^t, x) be a solution of the backward heat equation 

¥>t = -Vxx < t < T, X e R (2.3) 

with "final condition" 

Of course, v'(T, x) is the heat kernel evaluated at time T, and therefore, for 
< t < T, 

^{t. x) = I , exp --77-^ — - . 

^A'K{2T-t) V ^^T-t)J 

Next, a short calculation yields the following lemma. 
Lemma 2 Let ip{t,x) be as in ^7^). IfO<t<T, then 



{t,x)>V2ip^'^\T, 



ip^'^\t,x) ^ V'^f'^'[i,x 



Proof of Lemma |. Since < t < T, we have that T/(2T - t) > 1/2, and 
T-i > (2r -t)-^. Therefore 



^(^\t,x) Ft 



exp 



X 



2 



> V2 



T-1- {2T-t)-^ 



This proves Lemma §. 



For future use, we also compute the L^ norm of ipitjx)"", for a > 0. We 
claim that there exists a constant C = C{a) > 0, not depending on T, such 
that 



|y,(t,x)"||i = i^it.xYdx (2.4) 

.2 

= C'(a)(2T-t)(^^")/2 



(4.(27 -0)-"'^/exp-j^b. 



Now, let 

M{t) = I ip{t, x)u{t, x)dx, 0<t<T. 



and note that M{t) is a continuous J^t martingale for < t < T. This asser- 
tion was proven in Lemma 2.3 of [[MS93|] . One can also check it heuristically. 



by formally differentiating M{t) and applying ( |1 . 1| ) and ( |2.3| ). Lemma 2.3 of 
||MS93|| also states that < t < T, M{t) has square variation 



{M)t= / g{u{s,x))ip{s,xydxds > / I u{s,xY'^ip{s,xYdxds. (2.5) 

Of course, M{t) = M'^{t) implicitly depends on T. We now prove the fol- 
lowing lower bound on (M)^. 

Lemma 3 There exists a constant Ci > 0, not depending on T , such that if 
0<t<T, then 

Jo 

Proof of Lemma ^ 

Let 

Note that 

2-a , , 

^— + a = l (2.7) 

27 

and 

1 — a , ,1 , ^ 

__.(!_ 27) = --. (2.8) 



Furthermore, for t fixed, 



Mt,x)% 

is a probability density over x G R. Using Jensen's inequality, 
and (|2.8|) , we find that 

u{s,xY'^ip{s,xYdx (2.9) 

= Ms,xrh f u{s,xY^v{s,xr-^-^j^^^dx 

Jl \\Lp[s, xYWi 

> Ms,xrh (fuis,xMs,xf-y^'- \.^^''f.. dxT 
\Ji ||(/9(s,a;)"||i J 

= \\ip{s,x)"'\\i'' '^ i / u{s,x)(p{s,x)dx] 
= CiT-^/^M{sy\ 



where Ci = C{a)^~^'^, and a was defined in ( p.6|) . After integrating ( p.9| ) over 
s e [0, t] and putting this together with (|2.5| ), we get Lemma ||. I 

Using Lemma ^ it is possible to compare M{t) to a time-changed Brow- 
nian motion. In the standard way, the new time scale is given by {M)t. 
Let 

T(L) = 16^-2^8 

and consider the following gambler's ruin problem. Start with M(0) = 2. 
Let T = t{L) be the first time t that M{t) = 1 or M{t) = L. Using the 
optional sampling theorem in the usual way, we deduce that EM{t) = 2, 
and therefore (if M(0) = 2), 

P(Mir) = L) = j^. (2.10) 

In fact, we wish to show 

Lemma 4 //T = T(L) = IQC^^L^, then 

P{M{tAT) = L)> 



2(L 
8 



Proof of Lemma ^. 

The definition of r implies that for all t G [0,r], we have M{t) > 1. 
Therefore, by Lemma ^, if t G [0, r] then 

{M)t > CiT-^'H. 

Now M{t) is a continuous super martingale, so it follows that M{t) is greater 
than or equal to a time-changed Brownian motion with time scale {M)t. In 
other words, for some Brownian motion -B(t), we have M{t) > 2 + B{{M)t). 
Therefore, since 

{M)t > CT-^/'^T = CT^^\ 

we have 

P{T <t) = P (T < r < a, 1< M{t) < L for t G [0, T]) 

< P (T < r, 1< 2 + B{{M)t) < L for < t < T) 
= P (T < r, 1 < 2 + P(t) < L for < t < {M)t) 
= P (t < r, 1< 2 + P(t) < L for < t < CiT^/^ 

< P I sup P(t) < L - 2 

yte[o,CiTi/2] 

Using the reflection principle, we continue with 
P{T <t) 
< 1 - P I sup B{t) > L - 2 



1 - 2P [B{CiT^/^) >L-2 
P(|P(CiTi/')| <L-2 



Therefore, if 



< 2(L-2)(27rCiTi/2)-i/2 

T = T{L) = leC-^L^. 



then 



Then, by (|210|) , 



P{T <r) < 



1 



P{M{tAT) = L) > P{M{tAT) = L,T>t) 

= P {M{t) = L,T>t) 

= P{M{t) = L)-P{M{t) = L,T <t) 

> P{M{t) = L)-P{T <t) 
1 1 



> 



L-1 
1 
2(L- 



2(L 



This proves Lemma |. I 

Now we continue with the proof of Theorem |l|. Lemma § imphes that 



V 



> 



>(2L) 



-1 



2(L-1) 
Using Lemma ||, with L = L/2, we deduce that 

N > K-\L/2y^^-^^ - 1 > K~\L/Ay^^-^ 
if L is large enough. Thus we find that 



(2.11) 



if 7 > 3/2 and L is large enough. 

To finish the proof, we can apply the same argument as in | MS93 |, sections 
3 and 4. Since these arguments carry over, word for word, we will merely 
summarize the argument here, and refer the reader to [[MS93|| for details. 

First, we need to split up the solution u. For this, we quote Lemma 2.5 



of ||MS93|| . But first, define 

b{x,y) 



[x + yy-y - y 



27 



10 



Lemma 5 For t > 0, x G I, i = 1,2, ... ,N consider the N recursively 
defined equations 

u\t,Q) = u'{t,J) = 
u\0,-) = ul 

where u^ = by definition. Here the {W^} 's are independent white noises 
and the {uq} are some collection of nonnegative initial functions. Let us then 
define the process 

~ .^ ^ ^ f Ell u\t, ■) /or < t < min{o-(«0 ■ t = 1,2, . . . ,N] 
I cxD otherwise 

for all t > 0. Here, cr(-u*) denotes the blow-up time a with respect to m*. For 
t > 0, X G I, we have that u is a solution of 

du d'^u ^^^-^ 

— = h u'W 

dt dx^ 

M(t,0) = M(t, J)=0 

N 

for some white noise W which is a linear combination of the {W^}. 

We use Lemma ^ to split up the solution u into the sum of solutions m*. 
Later, we will further split up the u*. Section 4 of ||MS93|| explains how to 



use Lemma |^ to split up u over and over again, at a sequence of stopping 
times. Each of these smaller solutions will have a larger noise term than in 
(|1.1| ), so the corresponding total mass martingales U^{t) = JiU^(t,x)dx will 
have 

{U')t> t U\sf^ds. 
Jo 

We need a way to split up u, given that a certain integral is sufficiently large. 

The following lemma is an easy modification of Proposition 3.2 of |[MS93 |. 



Lemma 6 Let E'^{J) denote the class of nonnegative C°° functions on 
[0, J]. There exists a constant K > {) such that the following holds. Let 

11 



J > 4: be fixed. Set J = 72^(^-1) . If N > is and integer, and /o G E^{J) 
satisfies 

f^ cl> (t, x22(i-^); ^0, j) fo{x)dx > KN, 

for some zq in [1, J — 1] and some < t < 1, then there are functions 
{/i : z = 1, 2, . . . , A^} C E^{J) such that 

N 
i=l 

and for each i = 1,2, . . . ,N, 

rJ 







(t){Q,x]Zi,J)fi{x)dx>2. (2.13) 



for some z.i in [1, J — 1] 



In | |MS93| | , Lemma ^ was shown for A^ = [2^"^"^] , but the proof given there 
also imphes the above result. 

Now we continue the main argument. We can assume without loss of 
generality that 

j^G{T,x,y)u{Q,y)dy>2. 

If this condition fails, wait until time 1, when it has a positive probability of 
holding. Now wait until time T. By Lemma p^.lU| , we have that 

p(^JG{2T,x,y)u{Q,y)dy>L^ > -^^—- = p (2.14) 

Let 

N = K-^L^^'<~^\ (2.15) 

Now perform the scaling as in Lemma [^, with L = L/2. For the scaled 
function v, we see that 

f G{2T, X, y)u(fd, y)dy > L^^^-^^ = KN. 

Then, Lemma |^ shows that we can decompose 

N 

u{t,x) =Y,fiix) 



i=l 



12 



such that for some set of points {zi}f_-j^, 

'G{2T,x,y)My)dy>2. 



We use these fi as initial conditions for new functions M*(t,x), which satisfy 
(|2.12| ), and we call these M*(i, x) offspring of u{t, x). 



If 

' G{2T,x,y)u{0,y)dy<L, 



II 
then we say that mass has died. 

Repeating the argument, we find that there is mass alive at stage k if the 
branching process of the m's is alive at stage k. But this is a Galton- Watson 
process with expected number of offspring at least 

pN = K-'L^i-r-^^—^-— > 2-iir-iL2(^-3/2) 
2(L - 1) - 

by ( |2.14| ) and ( ^.151 ). Therefore, if 7 > 3/2 and L is large enough, the 
expected number of offspring is at least 

pA^ > 1 

and there is a positive probability of survival. But survival means that there 
is mass present at each stage. This, in turn, means that u{t, x) blows up in 
finite time. Therefore, there is a positive probability of finite time blow-up. 
But this conclusion contradicts our assumption (|2.1|) that P{a < cxd) = 0. 
Thus, Theorem |l] is proved. I 
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